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Introduction 

This is a request for renewal of a research project which is 
concerned with the development of new and improved methods for. 


• Multiresolution representation for very large, time- 
dependent, curvilinear grids with multiple masked 
overlapping blocks. 

• Explicit methods for streamlines, stream surfaces and other 
flow field visualization methods. 

In the following two sections we describe the current progress of 
the research project and a statement of the research work for next 

year. 


Report on Current Progress 

We first report on our current progress in the area of explicit 
methods for tangent curve computation. The basic idea of this 
method is to decompose the domain into a collection of triangles (or 
tetrahedra) and assume linear variation of the vector field over each 
cell With this assumption, the equations which define a tangent 
curve become a system of linear , constant coefficient ODE's which 
can be solved explicitly. There are five different representation of 
the solution depending on the eigenvalues of the Jacobian. The 
analysis of these five cases is somewhat similar to the phase plane 
analysis often associate with critical point classification within the 
context of topological methods, but it is not exacdy the same. There 
are some critical differences. Moving from one cell to the next as a 
tangent curve is tracked, requires the computation of the exit point 
which is an intersection of the solution of the constant coefficient ODE 
and the edge of a triangle. There are two possible approaches to this 
root computation problem. See Figure 1. We can expiess the tangent 
curve into parametric form and substitute into an implicit form for 
the edge or we can express the edge in parametric form and 
substitute in an implicit form of the tangent curve. Normally the 
solution of a system of ODE's is given in parametric form and so the 
first approach is the most accessible and straightforward. The 
second approach requires the "implicitization of these parametric 
curves. The implicitization of parametric curves can often be rather 



difficult, but in this case we have been successful and have been able 
to develop algorithms and subsequent computer programs for both 
approaches. We will give these details along with some comparisons 
in a forthcoming research paper on this topic. 


Parametric Curve 



Implicit Edge 


Ax(t) + By(t) + C = 0 
No Bracketing 

0 < t < ? 



Figure 1. Two approaches 


Parametric Edge 
F(x(t), y(t)) = 0 

Bracketing 

0<J<1 

to computing exit point. 


We now include some examples illustrating the use of our new 
methods. In Figure 2, we show the topological graph of a simple 
vector field with is specified by the two equations. Since we know 
the exact solution in this case, we can apply our methods and see if 
they work at all and also see how well they approximate the true 
solution. We can vary the size and resolution of the grid in order to 
asses the rate of approximation. In the top image of Figure 2 we 
have used a 5x5 Cartesian grid with diagonals yielding triangles and 
in the bottom image we have use a 15x15 grid. The example of 
Figure 3 shows that these methods can be used on any domain which 
can be triangulated! Here we have a two-dimensional flow ie 
defined over spherical domain which has been triangulated. In this 
particular case, the triangulations are obtained by subdivision and 
the right column is one more level of the subdivision than the left. 
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Each column shows three different rotations of the Earth. The blue 
curves represent the topological graph and the cyan curves are some 
additional tangent curves. In Figure 4, we show the topological 
graph of a two-dimension flow defined over a curvilinear grid which 
has an interior boundary. We include this example so that we might 
discuss the points of attachment and detachment on the interior 
boundarv. The computation of this type of critical point has not been 
rigorously treated in the past. Within this context these points are 
precisely defined and easily computed. Without these points, the 
linking algorithm will leave some unattached foci and nodes in the 
domain. Note the connection between a focus point and a point of 
detachment in the lower left region. The details of this type of 
analysis and the improved linking algorithm will be included in the 
research paper we have already mentioned. In Figure 5, we show 
some results for a three-dimensional vector field. This is really 
preliminary work and we need to do a lot more in this area so as to 
make the algorithms efficient and robust, but the figure does show 
promise and the possibility to extend these methods to higher 
dimensions. 
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U(x,v) = - 0.103 + 0.052x - 0.303y 

+ 0.038xy - 0.232x 2 + 0.612y 2 

V(x,y) = + 0.143 + 0.688x - 0.145y 

- 0.213xy - 1.030x 2 + 0.246y 2 
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We next report on our research progress in multiresolution 
modeling. Our work in this area can be broken down into two 
subareas: i) triangular grids with applications to arbitrary domains 
and ii) curvilinear grids. 

First we discuss curvilinear grids. In the original proposal we 
mentioned the idea of using multiresolution models for parametric 
curves in this context. The idea is based upon die observation that a 
curvilinear grid can be viewed as a parametric map and researchers 
have already developed techniques for using wavelets on parametric 
curves. See Figure 6. 



x 


Figure 6. Curvilinear grid as a parametric map 

Unfortunately there is a basic problem with this approach and that 
lies in the fact that not only do you obtain a lower resolution 
approximation to the flow, but you also get a lower resolution 
approximation to die boundary of die domain. This means that die 
domain will be "smoothed down" also. See Figure 7 and note diat the 
inside corner in die lower left region has been rounded off. 
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Figure 7. The boundary problem with the parametric wavelet 

approach 

After some consideration, we felt it was rather important to keep 
the domain intact because often in a flow analysis, the most 
interesting and important activity occurs near to and as a result of 
the boundary. In order to preserve tire boundary, we have adopted 
a domain decomposition scheme which is based upon alternate "knot” 
removal. This is illustrated in Figure 8. Up to now, we have 
developed and implemented some Haar type of wavelets based upon 
this type of domain decomposition, hi Figure 9 we show an example. 
This example has both an internal and external boundary and the 
reader will note that these boundaries stay intact throughout the 
different levels of decomposition. The example of Figure 10 
illustrates the capability to do region selective (lossless) 
reconstruction. The decomposition allows us to find the best lower 
resolution approximation and the reconstruction process allows us to 
add in "detail" when and where it is desired. The possibilities 
indicated by this rather crude example are really quite exciting! 
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Figure 8. Decomposition by knot removal 
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We now discuss our research progress on multiresolution methods 
for nested triangular grids. We have developed some new Haar type 
wavelets for subdivision of triangular domains as illustrated in 
Figure 11. The basic decomposition step is illustrated in Figure 12. 
The <j>'s are the characteristic functions for the subtriangles and 
therefore form a basis for piecewise con tan t functions. The y's are 
the wavelet functions which capture the detail (error) and hopefully 
they are mutually orthogonal so that we obtain best least squares 
approximations when we do the wavelet decompositions. This 
decomposition is simply a change of basis which can be represented 
by the refinement equations 


<5>(p) 


h 0 h a hb he 

Ya(p) 


gao gaa gab gac 

Vb(p) 


gbo gba gbb gbc 

- Tc(p) - 
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<t>a(p) 

4>b(p) 

- <j>e(p) J 


( 1 ) 
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Figure 12. Basic decomposition step 


As we have set up things here, we will have defined Haar wavelets 
for these types of subdivision domains once we specify the h s and 
g's of the refinement equation (1) above. If the aieas of the 
subtriangles T 0 , T a , T b and T c are all the same then it is possible to 
define Haar wavelets that are mutually orthogonal and affine 
invariant. If the areas are not the same, then we must settle for 
biorthogonal wavelets which means that we lose the best 
approximation properties. One applicadon area where we definitely 
have subtriangles with different areas is the spherical triangulation 
obtained by subdivision. Figure 13 illustrates this property. 


15 


oc 0 * a a ^ a b ^ a c 


Tetra, Level 1 




Level 3 ^ vel 6 . 

Figure 13. Spherical Triangular Grid 


In a 1995 SIGGRAPH paper, P. Schroeder and W. Sweldens 
presented some biorthogonal wavelets for spherical triangulations. 
The definition of their wavelets are given by the refinement 
equations 
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We have developed what we consider to be some improvements 
over these wavelets. These new wavelets are only biorthogonal also, 
but have the property that when the areas do become the same, the 
wavelet become mutually orthogonal. They are more nearly 
orthogonal and so we obtain better approximations with die wavelet 
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decomposition. The reader will note from Figure 13 that for 
spherical triangular grids these areas quickly become very similar 
and so for this application these new spherical wavelets will tend 
towards orthogonal wavelets. Ihe refinement equations for these 
new wavelets are 


<Kp) 

Va(P) 

Yb(p) 

-Yc(p)J 
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We have applied our new wavelets to data consisting of a vector 
field defined over the sphere. In the example of Figure 14, we took a 
known vector field function and evaluated it to obtain the data. In 
the left column we show three views of the reconstruction of the 
wavelet decomposition using only 1% of the coefficients. This means 
we started with a vector valued function which had 2048 coefficients 
(one for each triangle). The first time we apply the decomposition 
equations we obtain a function which is piecewise constant over 512 
triangles and plus the 1536 detail (or wavelet) functions. No 
information is lost, just a change of basis. The next step yields a 
lower resolution approximation consisting of 128 terms and o 84 
detail at this level plus the 1536 detail from the first decomposition 
step. This is continued two more times. This final expansion (of 
2048 terms) is scanned and the 20 (1%) functions with the largest 
coefficients are summed to obtain the partial reconstruction. In the 
right column we show the reconstruction based upon the largest 409 
(20%) coefficients. It is indistinguishable from the original function . 
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Figure 14. Partial tnearlv orthogonal) spherical wavelet 
reconstruction ot a (low lield de lined over a spheiical domain. 
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In this next example shown in Figure 15, we applied these new 
spherical wavelets to some "real" data. This is data we obtained from 
Roger Crawfis and Nelson Max of Lawrence Livermore National 
Laboratory. Actually it is simulated data resulting from a global 
weather model. This data represents one slice in elevation and one 
time step. The first column shows three different views using a 
topological graph of the original data. In the second column the 
reconstruction based upon the 3% largest coefficients is shown. This 
low resolution approximation allows us to get an overview of the 
flow without the clutter and detail. As it was shown with the 
example of Figure 10 it is possible with this type of analysis to zoom 
in and add detail to any region we wish. Most often wavelets and, in 
general, multiresolution models are used for compression or 
efficiency in fly-throughs. The use of these types of models for 
analysis and visualization of scientific data is potentially very useful 
and a very exciting area of research. 
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wealher model 
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Statement of Work for Next Year 


Concerning work in the area of explicit methods for computing 
tangent curves, we feel that our algorithms and computer programs 
are quite robust and efficient at this time. We want to publish tins 
research and make the code available with the hope of getting some 
feedback. We want to do some further work on comparing the new 
methods to the conventional stepping methods. As we have mention 
before these new methods allow errors to be bounded and controlled 
while all other methods only can at best estimate errors. We feel 
that this is important in general and in particular for linking 
algorithms for topological methods. We hope to determine just how 
critical this issue is. Even though we showed a 3D example in Figure 
5 this was really just preliminary work. We need to do a throug 
and complete job on the case analysis based upon eigenvalues. Also 
as as we pointed out in die 2D work, the approach which uses an 
implicit representation for the tangent curve and a parametric 
representation for the edge has some definite advantages. We need 
to do the mathematics of finding implicit representations of the 3 
parametric tangent curves. Once all the mathematics is done we 
need to develop and implement the algorithms. Linking algorithm s 
for topological methods require die definition of points of attachment 
and detachment. As we mentioned, in the past this has not been 
dealt with very well. For our linearly varying fields in 2D, can 
precisely define and compute tiiese special critical points. The same 
theory for 3D needs to be developed and we intend to do this. 


Concerning work in the area of multiresolution methods for 
curvilinear grids, we want to develop some methods which are 
higher order and have more continuity than the simple ^aar 
wavelets diat we have developed so far. We will first do this for 2D 
flows, starting with the development of piecewise linear wavelets 
over the nested decompositions we already mentioned. Also, we 
want to extend die basic 2D Haar wavelets to 3D. This requires the 
development of a 3D knot removal scheme leading to the nested, 
cellular decomposition of the domain. In principle this should not be 
much of a problem, but there are some details with the 
tetrahedrization of cells that we are concerned about but we feel that 
they can be dealt with. 


Concerning our work next year in the area of multiresolution 
models over triangular grids, we need to extend these basic Haar 
methods to higher order and higher continuity. Also we would like 
to extend this work to 3D. A critical component of the 2D 
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implementation was a clever triangle naming scheme which allowed 
neighboring triangles to be quickly and easily determined. To some 
this might be considered just an implementation detail, but to us^it 
was very critical to the entire algorithm. We need to develop a ^>D 
version of die labeling scheme and then proceed to develop die basic 
refinement and decomposition equations for the Haar wavelets 
defined over nested tetrahedral subdivisions. Just as the 2D 
triangular methods were extended to a spherical domain, we should 
be able to extend the 3D version to tetrahedral approximation of any 
three dimensional manifold. This would include SO(3) (die space of 
orientations) and could have some interesting implications on image 
interpolation and its use in VRML applications. 
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In conclusion, we give an outline of the planned research activities 
for next year: 


• Explicit Methods 

Publish papers, distribute algorithms, solicit feedback 
Extend to 3D 

Implicit representation of 3D parametric curves 
Definition of attachment and detachment curves on 
interior boundaries 

Tetrahedrization of 3D curvilinear grids 

• Multiresolution over Curvilinear Grids 

Two-dimensional 

Develop higher order (continuity) wavelets 
Apply techniques to many and various data sets 


Three dimensional 

Develop knot removal algorithm and software for 
3D curvilinear grids. Develop Haar wavelets 
over these grids. 

Develop higher order (continuity) wavelets 
• Multiresolution methods for nested triangular grids 


Two-dimensional 

Develop higher order (continuity) wavelets 
Apply to spherical data and flow over other 
domains 

Three dimensional 

Develop labeling scheme and basic Haar wavelets 
Apply to image interpolation problem for VRML 
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Budget 

We are requesting funds for the principal investigator as follows. 

Gregory M. Nielson 

2 summer months 

(see budget sheet for amounts) 

We are requesting funds for two students. 

Two ASU Graduate Students 

9 months half time, 3 months full time 
(see budget sheet for amounts) 

We are requesting funds in the operations category as explained on 

the budget sheet 
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ESTIMATED BUDGET 


PERIOD OF PERFORMANCE: January 1,1995 through December 31. 1997 


SALARIES & WAGES 

A. PI: Nielson 

B. Co-PI: 

C. Co-PI: 

D. Co-PI: 

E. Postdocs 

F. Technicians 

G. Graduate Associates 

H. Graduate Assistants 
]. Undergraduates 

J . Secretary 

TOTAL 



0.00 

pm/AY + 

2.00 

pm/ Sum 


0.00 

pm/ AY + 

0.00 

pm/Sum 


0.00 

pm/AY + 

0.00 

pm/Sum 


0.00 

pm/AY 4- 

0.00 

pm/Sum 

0 0 

0.00 

pm/C Y 



0 0 

0.00 

pm/C Y 



2 <2> 

9.00 

pm/AY + 

6.00 

pm/Sum 

0 0 

0.00 

pm/AY + 

0.00 

pm/Sum 

0 0 

0.00 

pm/AY + 

0.00 

pm/Sum 

0 0 

0.00 

pm/C Y 




Year 1 

Year 2 

Year 3 

Summary 

S 1 6,610 

S 17.607 

SI 8,663 

S52.SS0 

so 

so 

SO 

SO 

SO 

SO 

SO 

SO 

SO 

SO 

SO 

SO 

SO 

SO 

SO 

so 

so 

$0 

SO 

so 

$35 ,020 

S3 6, 050 

S37.767 

S10S.S37 

$0 

SO 

SO 

so 

SO 

SO 

SO 

so 

$0 

SO 

$0 

so 

$5 1 ,630 

S53.657 

S56.430 

$161,717 


FRINGE BENEFITS 

A. Faculty at 

B. Staff at 

C. Students at 

TOTAL 


S4.153 

so 
S 1,051 

54,402 
SO 
SI ,0S2 

$4,666 
SO 
SI, 133 

$13,221 

SO 

S3, 266 

S5.204 

S5,4S4 

S5,799 

S16.4S7 

SO 

SO 

SO 

SC 


III. TRAVEL 


[V. OPERATIONS 

A. Consultants 

B. Publication/Pace Charges 

C. Photocopy, telephone, postage, misc. 

D. Material & Lab Supplies/Equipment Use Fees 

E. Computer Software 

F. Lab Equipment under $500 

G. Subcontracts 

H. Participant Support Costs 

1 Tuition 0 © S1.S95 per AY 

J. Student Stipends 

K. Other 

TOTAL 

V. CAPITAL EQUIPMENT 

VI. TOTAL DIRECT COSTS 
V!l. INDIRECT COSTS 

Year 1 
Year 7 
Year 3 

Viil. TOTAL PROJECT COSTS 


(4 9c escalation) 


5: OF M 1 DC 
5 3 5 E M T D C 


SO 
SO 
S500 
SI ,000 
S200 
SO 
SO 


SO 

so 

S500 

SI ,000 

so 


$0 

so 

$500 

si.ooo 

so 


SO 

SO 

SO 

SO 

SO 

SO 

SI. 700 

SI .500 

SI .500 

SO 

SO 

SO 

S5S.534 

$60,64 1 

$63,729 

S30A3S 

S 3 1 . S 3 7 

$33. 45S 

SSS.972 

S92.47S 

$97. 1 $7 



